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Abstract. We have studied the atomic density of a cloud confined in an isotropic harmonic trap at the 
vicinity of the Bose-Einstein transition temperature. We show that, for a non-interacting gas and near this 
temperature, the ground-state density has the same order of magnitude as the excited states density at the 
centre of the trap. This holds in a range of temperatures where the ground-state population is negligible 
compared to the total atom number. We compare the exact calculations, available in a harmonic trap, to 
semi-classical approximations. We show that these latter should include the ground-state contribution to 
be accurate. 
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The phenomenon of Bose-Einstein condensation (BEC) 
is a phase transition. Below the critical temperature T c , 
the ground-state population, which is the order param- 
eter, becomes macroscopic. This phenomenon, that hap- 
pens strictly speaking only at the thermodynamic limit, is 
usually illustrated in textbooks with a homogeneous gas. 
Experimentally, the Bose-Einstein condensation of dilute 
gases has been observed since 1995 with atoms confined in 
a harmonic trap . These stimulating experimental data 
have quickly pointed out that two effects had to be taken 
into account: the interatomic interactions and the finite 
number of atoms 0- Several papers, as the present one, 
have studied harmonically trapped ideal gases containing 
a finite number of atoms. Two quantities have been in- 
vestigated in detail: the atom number |3l4lo17l8l9"] and 
the specific heat |5I7I9) . For a finite but large (typically 
10 6 ) number of atoms, the properties of the atomic cloud 
change abruptly at a characteristic temperature we will 
name the transition temperature T* . This temperature is 
shifted compared to T c , but by a small amount, typically 
of few percent for atom numbers around 10 6 . There is also 
a characteristic temperature for the specific heat; it is dif- 
ferent from the previous one but still close to T c 0*3]. 

Surprisingly, less attention has been paid on the atomic 
density of an ideal gas 1(} . In a homogeneous gas it is ob- 
viously equivalent to the atom number but this is no more 
the case in a spatially varying potential. It becomes the 
good parameter of the theory, in particular to perform 
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local density approximations. This quantity is then par- 
ticularly important for the study of the shift of the critical 
temperature by the interatomic interactions, both within 
the mean-field approximation and beyond this approx- 
imation |TJ|. We will show, in the case of an isotropic 
harmonic trapping and for a finite atom number, that the 
ground-state density at the centre of the trap increases 
much more sharply than its population as the tempera- 
ture decreases. This leads to the fact that near the Bose- 
Einstein transition temperature the density is already dom- 
inated by the ground-state contribution. This holds what- 
ever the atom number is, and is a remanence of the infinite 
compressibility of an ideal gas at the thermodynamic limit 
|12j . Usual semi-classical approximations do not take into 
account the ground-state contribution and then fail in the 
vicinity of the Bose-Einstein transition temperature. This 
is not a finite size effect in the sense that it is not related 
to the discretization of the excited states energy levels. We 
will compare the exact results with semi-classical approx- 
imations. The addition of the ground-state contribution 
on the latter ones improves their accuracy. We will finally 
show that the influence of the ground-state is smaller if 
the measured quantity is the density integrated over at 
least one dimension. It is still large for typical experimen- 
tal parameters. 

We will perform our calculations in the grand canon- 
ical ensemble (GCE). Then, the Bose-Einstein distribu- 
tion gives the population iV, of a given energy level ef. 

Ni = (e 13 ^-^ - l)- 1 with N i = N - Here P = 1 / k BT 

i=0 
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with ks the Boltzmann's constant, p the chemical poten- 
tial and iV the total atom number. The equivalence be- 
tween GCE and the canonical or microcanonical ensem- 
ble, these latter being probably more appropriate descrip- 
tions, is generally not guaranteed, especially for systems 
that are not at the thermodynamic limit. For instance, it 
is well known that the GCE predicts unphysical large fluc- 
tuations of the condensate population at low temperature 
[15]. However, the authors of Ref. [10114115] have shown 
that the occupation numbers Ni in GCE are very close 
to the ones in the canonical ensemble. The difference is 
more pronounced for small atom number and anisotropic 
clouds. As a result and because GCE enables to give an- 
alytic expressions on contrary to the other ensembles, we 
will use GCE in the following. 

For a fixed atom number, the chemical potential in- 
creases as the temperature decreases. As p has to be smaller 
than eo, the ground-state energy, the excited states pop- 
ulation will saturate when p approaches eo whereas iVo is 

oo oo 

still increasing: N - N Q = £ Ni(p, T) < £ JV,-(e , T). As 

i=i i=i 

in Ref. |2I16| , we will define the transition temperature T* 
as the temperature for which the excited states saturated 
population is equal to the total atom number: 

oo 

Y / N i (e ,T') = N (1) 

i=l 

As pointed out in the introduction, there is not a unique 
definition of the transition temperature for a finite atom 
number. Other definitions use, for instance, a change in 
the slope for the condensate fraction in function of temper- 
ature (more explicitly d = 0) [TJ\, a change in the 
power dependence on the condensate fraction in function 
of the atom number which are also pertinent. We have 
checked that these various definitions affect marginally the 
value of T* and do not modify our conclusions JS| ■ I n the 
following we will then use Eq.QJ to define T*. Note that 
the chemical potential p* at the transition temperature 
is close but not equal to the ground-state energy; it is 
determined by the constraint 

00 

J2Ni(»*,T*)=N (2) 
t=o 

There are only a few examples of trapping potentials 
where the eigen-energies and the eigen-functions are known 
exactly. Semi-classical approximations give usually accu- 
rate enough results and are suited to include interatomic 
interactions, at least perturbatively. We will derive vari- 
ous type of semi-classical approximations in the following 
and test their accuracy because the harmonic potential is 
an exactly solvable potential. 

We will first examine the situation where huj < fcgT 
with lu the oscillation frequency of the isotropic harmonic 
trap. This corresponds to the large atom number limit 
and semi-classical approximations should work. Replac- 
ing the discrete energy spectrum by a continuous one 
and neglecting the ground-state energy the density is 



p(r) = j3-(7| [z exp(— ^(r/a) 2 ] with z — the fugacity, 
t = -j^Hp and g& () a Bose function J^I- With the above 
notation, the thermal de Broglie wavelength is A = o^JI-kt 
and the size of the cloud Similarly, the 

atom number is N = g^(z)/T 3 . Equation (JIJ leads then 
to N = C(3)/r* 3 , with t* the value of r at T = T* . The 
above expressions for the density and atom number are 
in fact approximations for the excited states and do not 
contain the ground-state contribution. Then p* defined by 
Eq.(J2J is equal to and z* = 1. The transition tempera- 
ture defined here corresponds to the critical temperature 
T c . The peak density at the transition temperature is then 
given by p(0)X 3 = gz{z*) = ((3/2) » 2.612. For temper- 
atures below T c , the excited states population is given 
by £(3)/t 3 . Then, the ground-state population fraction is 
N /N = for T > T c and N /N = 1- (T/T c ) 3 for T <T C . 
This fraction will be plotted in fig^ labelled with sc^ . 

These approximations are too crude and give inaccu- 
rate results for the atomic density, however. The reason is 
that the ground-state contribution cannot be neglected. A 
better expression is p{r) = -p-ffa [ze^^^/ 17 ^ ] + po(r) and 

similarly N = ^g 3 (z) + N Q with p (r) = j^e^l^ 
and No = The value of T* is unchanged as it is 

defined by the excited states saturation, but z* is now 
different from 1. Using 53(2*) « C(3) — C(2)^* with z* — 
e~ x ' (x = (3(e — p) > 0), one finds using Eq.@ that 

x* « r* 3/2 /v/C(2) M The ground-state population is 
~ 1/x* and, as expected, is vanishingly small as r* — * 
compared to the excited-state population C(3)/r* 3 . The 
ground-state peak density is ~ jTT^pjs^i whereas the ex- 
cited state peak density is ((3/2) /A* 3 . As A* = u\/2-kt* , 
the two quantities have the same order of magnitude! 
The above high-N analysis predicts then that the degener- 
acy parameter at the transition temperature is p(0)A 3 = 
C(3/2) + 2^2C(2) w 6.24 and not 2.612. The ground-state 
population is extremely small but the size of its wave- 
function is also extremely small compared to the atomic 
cloud size. For a harmonic trap both depend on the same 
small parameter, raised to the same power. So, even for 
very large atom number, the traditional criterion for BEC 
should be modified. This effect is linked to the patholog- 
ical behaviour of the ground-state density at the thermo- 
dynamic limit, i.e. the infinite compressibility of an ideal 
gas This limit means N — > 00 with Nut 3 — > constant. 
The ground-state size being a — y/h/muj, the density of 
that state behaves as below threshold and is then 
infinite at the thermodynamic limit whereas the density 
above T c is finite. 

We will now address the case of atom numbers in the 
accessible experimental range, 10 3 — 10 6 . It is well known 
that the transition temperature will be shifted compared 
to T c 1I7| . A better approximation, which takes into 
account the ground-state energy to first order, is p{r) = 
-^{ 5 3[z(r)] + %gi[z(r)]} where f(r) = zeS (r/<T)2 . Then 
N — -\[g3(z) + 4r<?2 (z)]. The corresponding transition 
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temperature is T* c such that N = -MC(3) + §C(2)t* c ]. 
This is the usual semi-classical approximation found in the 
literature. The ground-state population fraction is then 
N /N = for T > T* c and N /N = l-(^) 3 «3)+^C(2) 

for T < T* c . This fraction, also plotted in ng(I] wm be la- 
belled with scq. Note that gi(z) diverges at z = 1 |2U) . 
meaning that this approximation is intrinsically inaccu- 
rate near the centre of the trap and near the transition 
temperature. This divergence is however weak, and any 
spatial integration would give a finite result. We can still 
cure this pathology by adding, as before, the ground-state 
contribution. We obtain then 



P S c(r) = ^{ 5i [5(r)] + f<?i[z(r)]} + ^ 



T* such that N = 



T92(Z)] + 

= tMC(3) 



l-z 
_i_ 3 



C(2)r* 



(3) 

This semi-classical approximation will be labelled with 
sc in the following. The comparison of T* c with the value 
given by the exact model (see below) can be used to check 
the finite size correction. Even so, this comparison is use- 
less to check the contribution coming from the ground 
state since it does not depend on it (same transition tem- 
perature as scq). 

We can now test these semi-classical approximations 
for a harmonically trapped gas. As we referred before, 
for this case, the eigen-energies and the eigen-functions 
are known exactly. The corresponding expressions of the 
atomic density and atom number labelled with ex in 
the following, are : 




"i 1 r 
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Fig. 1. Ground-state population fraction in function of the 
temperature in Tuo/kB unit for a cloud of 10 3 atoms. The dot- 
ted curve corresponds to the exact result given by model ex. 
The solid, dot-dashed and dashed lines correspond respectively 
to the semi-classical models sc, sco and SCoo- The last two ne- 
glect the ground-state contribution above their corresponding 
transition temperature and the first two take into account finite 
size effects. The model sc is the closest to ex near Bose-Einstein 
transition. 
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Pex(r) 



(v /^)3 A. (1 



r!)3/2 



i=i 



T* ex such that TV = E(7T-^ 
(=1 K ' 



-tanh(^)(J) 2 



1) 



where, here z = e 13 ^ 1 e °K The semi-classical model corre- 
sponds to a Taylor expansion in r of these last expressions. 



Fig. 2. Relative shift of the semi-classical transition temper- 
atures T c (dashed line) and T s * c (dotted line) to T£ x (see text) 
in function of the atom number. Both temperatures converge 
for high atom numbers. The critical temperature at thermody- 
namic limit, T c , deviates by less than 1% for N > 5 f . The 
semi-classical transition temperature denned for a finite atom 
number, T* c , is much more accurate and deviates by less than 
1% for N > 400. 



In fig^ we plot the ground-state population fraction 
in function of the temperature for the various models de- 
scribed above. When the number of atoms is only 10 3 , 
finite size effects are large. The prediction of model sc^ 
is clearly wrong compared to the exact model prediction. 
On contrary models scq and sc give a result close to the 
one of ex [21] • Figure [3 shows the relative deviations of 
T c and T* c from T* x in function of the atom number. As 
expected the different values are similar but, as above, the 
model sc give a closer result to ex than SCoo . The value T* c 
deviates less than 1% for N > 400 and the relative shift 
is ~ f0~ 4 for typical experimental atom numbers. This 
is well below actual experimental uncertainties. The ther- 
modynamic value T c deviates more, typically 1 % but is 
still close to T* x |3I4I7I9| . The discrepancy with T c would 
have been more pronounced for an anisotropic trap (see 
below) . 



This two figures illustrate what is called finite size ef- 
fects, the fact that the energy level spacing is not negligi- 
ble compared to the temperature. What we are interested 
in is the role of the ground-state. For this, the transition 
temperature and the condensate population fraction are 
not the best observables. It is nevertheless already clear 
from fign that sc is a significant improved model to de- 
scribe semi-classically a cloud near degeneracy compared 
to scq. The high-N model predicts that the ground-state 
influence should be much more pronounced on the peak 
density. We will now focus our attention on that observ- 
able, only in the more pertinent comparison between the 
models sc and ex. 

This is first illustrated on fig. [3] where the degeneracy 
parameter p(0)A 3 is plotted in function of the atom num- 
ber for clouds at T = T*. We plot this number for the 
semi-classical approximation sc and for the exact model, 
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Fig. 3. Degeneracy parameter p(0)A 3 in function of the atom 
number N for clouds at the transition temperature. The dot- 
ted line corresponds to the semi-classical model sc at T = T s * c 
and the solid line to model ex at T = T* x . Even if the degener- 
acy parameters are somewhat different, they both differ signifi- 
cantly to the usual value of 2.612 (dashed horizontal line). This 
deviation is due to an under-estimation of the ground-state 
density. The actual values are close to our high-N prediction 
of 6.24 (see text). 



ex. The two curves are higher than 2.612. This highlights 
the inaccuracy of the standard semi-classical models (sco 
or sCqo) that do not take into account the ground-state 
contribution. It confirms also the calculation developed 
above. The degeneracy parameter is astonishingly con- 
stant till 10 3 atoms and does not differ much even for 
smaller atom numbers. Models sc and ex, which have 
almost the same transition temperature, have the same 
asymptotic value of the degeneracy parameter. This value, 
6.24, is the one predicted by our high-N analysis. The 
model sc is significantly higher than this value for exper- 
imentally accessible atom numbers. This is because our 
first analysis does not take into account the §r term of 



model sc. To first order JHj , x* s 
and is then slightly smaller than 




. Consequently the 



ground-state peak density is bigger at T* c using model sc 
than at T c using the high-N model. The excited states 
peak density is also higher in model sc because of this 
§T<?i term. 

The next three figures deal with the cloud properties 
around the Bose-Einstein threshold. Figure 0] and figOD 
show the evolution of the condensate fraction N /N and 
the condensate peak density fraction in function of T for 
two different atom numbers, 10 6 and 10 3 . Figure shows 
the density profile of clouds near degeneracy. What pre- 
vails in fig0]is the sharp increase of the condensate peak 
density compared to the condensate population. Moreover 
the models sc and ex give very close results validating 
our analysis on the ground-state contribution near degen- 
eracy. This means that the peak density is a much bet- 
ter marker of the Bose-Einstein threshold than the atom 
number. This feature is in fact used experimentally: the 
appearance of a small peak over a broad distribution is 
the usual criterion to distinguish clouds above or below 
the transition temperature. This sharpness also explains 





1.0 -i 


c 




o 


0.8- 


1 








' 

o 


0.6- 


to 




w 
c 


0.4- 


CO 




"D 




c 

o 


0.2- 


o 






00- 



— 1 — I — I — I — I — I — I — I — 1 

88 90 92 94 96 
Temperature (h.o. unit) 



100 



Fig. 4. Condensate atom number fraction No /N (dashed line) 
and peak density fraction po(0)/p(0) (solid line) in function 
of the temperature in harmonic oscillator unit hm/kB, using 
model ex. The cloud contains 10 6 atoms. The transition tem- 
perature is T* x — 93.37hui/kB and the asymptotic thermo- 
dynamic temperature is T c = 94.057iaj/fcs. The positions of 
these temperatures are shown as vertical lines in the figure. 
The ground-state peak density increases much more sharply 
than the ground-state population around the transition tem- 
perature. The former has also a significant value above T£ x . 
The model sc is indistinguishable for No/N, but is slightly 
different for po(0)/p(0) (dotted line). 
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Fig. 5. Same as in Fig^]but with 10 3 atoms. The transition 
temperature is T* x — 8.71hui/kB and the asymptotic thermo- 
dynamic temperature is T c = 9.41hui jhs- Since the number 
of populated states is considerably reduced compared to fig0| 
the discrepancy between sc (dotted lines) and ex is more pro- 
nounced. This also explains why the increase of the condensate 
peak density is slower. 



why the value of the peak density is very sensitive to the 
value of the temperature (cf. figEJ ■ Figure 0] shows also 
that, above threshold, the ground-state peak density frac- 
tion decays slowly. This is even more pronounced in fig El 
where N = 10 3 instead of 10 6 . It comes from the fact that 
the number of populated states is not macroscopic any- 
more (ksT < IOTiuj) and then the transition is smoother 
for smaller atom number. Once again, the density is a 
better marker of degeneracy than the atom number. This 
figure shows also that the |r term and the ground-state 
contribution make the model sc still very close to model 
ex, respecting the density and population fractions, even 
for 10 3 atoms. 

The above analysis is focused on the peak density i. e. 
at the centre of the cloud. Figure shows the total density 
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Fig. 6. Atomic density pc X in function of rjo where a is 
the size of the harmonic oscillator ground-state. The temper- 
ature is T = 93.37huj/kB and the atom number N spans from 
0.990 10 6 to 1.004 10 6 by step of 2000 atoms. The curve at 
threshold is in dotted line and corresponds to 10 6 atoms. The 
inset shows the excited states and ground state density profile 
at threshold. The dip around r = is mainly due to the first 
excited state population. 



profile of clouds, all at the same temperature, but contain- 
ing different numbers of atoms around N* x , the number 
of atoms for which T = T* x (N = N* x corresponds to the 
dotted line) . This figure simulates somehow an experimen- 
tal observation of BEC threshold. Only the central part 
is sensitive to the atom number; this corresponds to the 
condensate growing as the number of atoms is increased 
and to the fact that the excited states are already satu- 
rated for these atom numbers. Moreover, by looking at the 
graph, one would rather think that the Bose-Einstein tran- 
sition occurs for a smaller atom number. This points out 
that the definition on the transition temperature based 
on an atom number criterion does not fully correspond 
to the one based on the atomic density which would be 
more connected to experiments. The inset shows the ex- 
cited states and ground state density profiles at threshold. 
The excited states density exhibits a dip in the centre of 
the cloud, obviously not present in semi-classical models 
(monotonic functions). We check that the height of the 
dip is proportional to 1/t and can almost be totally at- 
tributed to the first excited state population. The aim of 
this paper is to show the importance of the ground-state 
in the study of non-interacting clouds close to threshold. 
The inset reveals that the first excited state density is 
also largely under- estimated; it represents ~ 10 % of the 
peak density whereas it contributes only to ~ 0.1 % of the 
population. 

We have shown results on the atomic density at the 
vicinity of the transition temperature. Detection techniques 
consist rather on ID-integrated density, corresponding to 
2D absorption images, or 2D-integrated density [53]. One 
can show that, at threshold, the ID and 2D-integrated 
peak density of the ground-state are vanishingly small 
for large atom numbers on contrary to the non-integrated 
case. The peak ID-integrated density fraction behaves at 
threshold as y^r and the 2D-integrated peak density as r. 
For typical atom number this is nevertheless not negligi- 
ble. This is illustrated in FigH where is plotted the con- 
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Fig. 7. Contribution of the ground-state on the peak density 
for, from bottom to top, ID, 2D and 3D images in function of 
the number of trapped atoms. The clouds are at the transition 
temperature T* x and the calculations use model ex. A 3D im- 
age would give the density in all three dimensions of space |25| 
whereas a 2D (resp. ID) image corresponds to the density in- 
tegrated over one (resp. 2) dimension. For = 10 4 atoms the 
ground-state contributes to ~ 26% in 2D images and ~ 6% in 
ID images. In contrast to 3D image, the ground-state contri- 
bution is very small for large atom number; it is not for typical 
atom numbers accessible in experiments. 

densate peak density fraction for 3D, 2D and ID images 
of clouds at threshold. The calculations use the model ex. 
At the transition temperature T* x , the ground-state con- 
tributes to more than 10% for N < 2500 atoms in ID 
images and for N < 8 10 6 atoms for 2D images. It means 
that, even with the conventional technique of absorption 
images, the effect should be experimentally observable if 
interactions could be switched off using, for instance, the 
magnetic tunability of the scattering length close to a Fes- 
hbach resonance |24j . 

Apart from the atomic density, two- and three-body 
inelastic loss rates will also be affected and could be 20 
to 30 % higher than predicted by model sen around the 
transition temperature for typical atom numbers. Finally, 
in most experimental set-ups, the trapping potential is 
anisotropic and finite size effects are then stronger. In- 
deed the term |r in Eq.© should be replaced by §=t, 
with u) — (Il^i) 1 ^ 3 the geometric mean and ui — ^ J^Wj 

i i 

the arithmetic mean j3j- Whatever the anisotropy is, & is 
always larger than lo, making the finite size contribution 
stronger. To first order and if kBT* x 3> hu>i for i = x, y and 
z, the ground-state contribution should be the same since 
our high-N analysis does not depend on any anisotropy. 

In conclusion, we have shown that the density of an 
ideal atomic gas is dominated by the ground-state contri- 
bution near the transition temperature. The inter-atomic 
interactions have been neglected in our analysis and will 
modify our conclusions. With repulsive interactions, the 
clouds tends to decrease its density at the centre of the 
cloud whereas it tends to increase it with attractive in- 
teractions. Previous calculations have treated separately 
finite size and interactions effects, both corrections be- 
ing finally added 0- Since the ground-state has a non- 
perturbative effect on the density, our analysis tends to 
prove that both effects have to be investigated together. 
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The approach of Ref.JIjOJ could in this respect provide 
helpful informations. Feshbach resonances, which enable 
to tune the interactions strength, constitute a powerful 
tool to check the accuracy of the different theoretical mod- 
els. Moreover, a full three-dimensional density measure- 
ment would also be valuable; this type of measurement is 
at the edge to be available in our experiment on metastable 
helium in Orsay 25 . 

We thank S. Giorgini for stimulating discussions. The Atom 
Optics group of LCFIO is member of the Institut Francilien de 
Recherche sur les Atomes Froids (IFRAF). 
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